Room temperature Rydberg Single Photon Source 
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We present an optimal protocol to implement a room temperature Rydberg single photon source 
within an experimental setup based on micro cells filled with thermal vapor. The optimization of a 
pulsed four wave mixing scheme allows to double the effective Rydberg blockade radius as compared 
to a simple Gaussian pulse scheme, releasing some of the constrains on the geometry of the micro 
cells. The performance of the optimized protocol is improved by about 70% with respect to the 
standard protocol. 
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The quest to find the perfect hardware to support 
quantum information processing is one of the main open 
problems to be solved to bridge the distance between 
quantum information theory and technological applica- 
tions. Indeed the conditions to be met are very de- 
manding and can be summarized in the need of having 
a (many-body) quantum system perfectly isolated from 
the environment but under perfect control at will. One 
of the promising candidates to fulfill such stringent con- 
ditions arc Rydberg atoms as their interaction strength 
can be tuned by twelve orders of magnitude by simply 
exciting them This very peculiar property is very 
difficult to be found in other architectures, as interac- 
tion is usually either very weak or always-on, and al- 
lowed for a fast development of their application in quan- 
tum information: phase gates 0, Q and recently a single 
photon source based on Rydberg excitation blockade [4| 
have been demonstrated using ultracold atoms. Despite 
these exciting successes, the need for ultracold tempera- 
tures still limits their applicability and thus recent efforts 
have been made towards replacing ultracold atoms with 
thermal vapor Q. Indeed there has been considerable 
progress toward similar results, demonstrating appropri- 
ate coherence times (gL four wave mixing [7J and van- 
der- Waals interaction [8j . Thermal vapor cells have been 
shown to act as a long lived quantum memory 0] as well 
as a probabilistic single photon source Thermal va- 
por cells have even been shown to support entanglement 
between quantum states of the spin variables in two dis- 
tant cells [10| . A global approach that combines quan- 
tum memories, single photon sources, quantum logic op- 
erations and efficient detectors on one platform is highly 
desirable e.g. to implement quantum repeater protocols 
like the one proposel by DLCZ (llj . 

Here we propose a scheme to obtain a single photon 
source using a thermal vapor of Rubidium 87 and we 
study theoretically its feasibility. The proposal stems 
from previous works originally thought for ultracold 
atoms where the Rydberg blockade has been exploited 
to prepare a collective single excitation state (a W-state) 
as a resource to accomplish different tasks [l^, EH ■ When 
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FIG. 1. Level Scheme of 87 Rb and sketched interaction. The 
used states are \g) — 5 2 Si/2, \i) = 5 2 Pa/2, \ r ) — 44D, |e) = 
5 2 Pi/2- Interaction occurs between Rydberg (Jr)) states with 
potential C 6 /r 6 and C 6 = h ■ 5880 MHz /im 6 [2<|. 



the W-state decays it will emit a photon with direction- 
ality given by the imprinted laser p hases. This decay 
process has been widely studied by 



14H18I] leading to 



the experimental realization of an ultracold single pho- 
ton source in 0]. We exploit optimal control theory to 
achieve the best possible preparation of the W-state un- 
der the experimental constraints in the more challenging 
regime at room temperature. Quantum optimal control 
theory has been applied successfully to solve this class of 
problems in few-body quantum systems in many different 
setups [Til and recently it has been proven to be success- 
ful also to optimally drive many-body quantum systems 
dynamics [lo, HH. We exploit this recently introduced 
technique to optimize the state preparation, simulating 
the unitary dynamics of the atomic ensemble driven by 
two lasers by means of the time-dependent Density Ma- 
trix Rcnormalization Group in its Matrix Product State 
(MPS) formulation that allows to include long-range in- 
teractions [22M24J I . This is the first application of opti- 



mal control theory to many-body quantum systems with 
long range interactions. We start our analyses with an 
ideal scenario of a frozen ordered unidimensional sys- 
tem: under such conditions, optimal control theory al- 
lows to find the best pulse to reach the W-state with 
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very high fidelity. We then extract the general features 
of the found optimal pulses using them as initial guesses 
to drive random instances of atoms distributed within 
a three-dimensional sphere and optimize them. The re- 
sulting pulses provide a robust strategy to prepare the 
desired state despite the new disordered geometry and 
system dimensionality. Finally, we study in details the 
effects of finite temperature (e.g. Doppler shifts) on the 
properties of the single-photon source and we show that 
we obtain a room temperature single photon source with 
performances similar to those obtained in the cold atoms 
case. 

System setup - We consider a hot rubidium vapor in 
a wedge shaped micro cell as recently introduced in 0, 0] . 
At room temperature we can expect a linewidth of about 
2 MHz for the 44D-state @. Fig. [1] shows the rele- 
vant levels of the 87 Rb atoms, where fli are the Rabi 
frequencies of the lasers and 5 g i, Si r and S re are the de- 
tunings of the corresponding transitions. By detuning 
f^i from the \g) <-> \i) transition we can adiabatically 
eliminate the level \i) in our theoretical model assuming 



. From now on we consider the three 
and |e) and two driving lasers with 
- QlQ z and n re = n 3 . The den- 
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level system \g), \r 

Rabi frequencies Cl gr = — 2 ~ s - anu a re = 0,3 
sity of the vapor can be varied by changing the gas tem- 
perature. At 220° C this gives a density of 87 Rb (27.83 % 
of all Rubidium atoms) of 543 atoms///m 3 . The veloc- 
ity distribution is Gaussian, corresponding to a thermal 
distribution. The decay rates are t\ = 26.2 ns for the 
\i) O \g) transition, and r 2 = 27.7ns for the |e) O \g) 
transition [28j | . As the single photon shall be emitted by 
the latter transition, this sets also the time scale for the 
decay process. This also limits the duration of the excit- 
ing laser pulses, that we set to be To = 2.5 ns. Given the 
separated timcscales with respect to the observed coher- 
ence times of about a hundred ns, we simulate a closed, 
coherent system as far as the state preparation is con- 
cerned and analyze the non-coherent part of the system 
evolution separately. Finally, the Hamiltonian of the sys- 
tem H = JT Hioc,i + Hi n t after the adiabatic elimination 
and choosing tti — f2 2 is given by 
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If we change the pulse duration from To to T and at the 
same time we scale the Rabi frequencies and detunings 
by Tq/T, while typical length scales like the diameter Lq 
of the system scale by Lt = Lq(T '/To) 1 / 6 the Hamilton 
practically remains the same. 

State preparation - The goal of the manipulation pro- 
tocol is to prepare a collective single excitation state 
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FIG. 2. Final infidelity e for a Gaussian 7r-pulses (circles) 
and optimal pulses (squares) in a one-dimensional N atoms 
chain with lattice constant ao = 0.35 /im at pulse duration 
To = 2.5ns and total extension Lq = ao(N — 1). 



means atom i in state e and all other atoms in state 
g), starting from the system's ground state \if>(0)) = 
\gg . . . gg) . In order to engineer the preparation of the 
W-state we simulate the time evolution of the system 
and optimize the pulse sequence. The time evolution is 
obtained by a tensor-network al gori thm including two- 
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body long-range interactions [22 
is performed by means of the Chopped RAndom Bases 
(CRAB) optimal control technique for many-body quan 
turn systems dynamics recently introduced in [ 
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The method is based on the a priori reduction of the 
problem complexity via a proper truncation of the acces- 
sible space used to describe the system wave function and 
the control field. In particular, in the system introduced 
above, the Rabi frequencies of two driving lasers fl gr (t) 
and Q er (t) are subject to optimization by the CRAB op- 
timization method. The pulses are expanded in a trun- 
cated basis fl; = c ijfij(t) wnere t ne c,j are the co- 
efficients of the expansion and the basis functions fij(t) 
are chosen according to physical properties of the system. 
In particular here we choose the first principal harmonics 
of the Fourier expansion. The figure of merit to be min- 
imized is the final infidelity of the system state evolved 
under the pulses O with respect to the target state \W), 
namely 



= l-|(IFhKT)}| 2 
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is 



\W) 



1 ^ w 



(where |e,-) = \g...geg...g) 



The multi- variable function defined by e = e(c. 
then minimized by means of direct-search methods 
Typical parameters are Ni = 14, 10 4 iterations with 10 4 
Trotter steps and bound dimension 10 : 40 for the MPS 
simulation. 

Results - We first concentrate on a one-dimensional 
chain of N Rydberg atoms with fixed positions. Despite 
its simplicity - as we shall show later - this model already 
captures the essential features of the system dynamics as 
the long range interactions effectively map this system 
in a particular instance of the full 3D scenario that we 
consider later. An experimental realization can be with 
optical lattices like in [2^. The interaction between two 
excited Rydberg atoms can act as a blockade for double 
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FIG. 3. Left: Gaussian 7r-pulses guess. Middle: Pulse from 
optimization of 4 atoms chain. Right: Guess optimized for 
3D cloud. Red: Q. gr , blue: Sl re . 



excitations and it enhances the driving Rabi frequency for 
the \g) <H> \r) transition by a factor of y/~N in an ensem- 
ble of N atoms, as one excitation suppresses N — 1 other 
possible excitations HE S3]. Fig. [Ishows the final infi- 
delity e with respect to the \ W) state obtained via Gaus- 
sian 7r-pulses (Fig. |3l left) of atoms in a ID chain with 
lattice spacing ao = 0.35 /im for To = 2.5 ns and total ex- 
tension Lq = ao(N—l). Notice that for different numbers 
of atoms TV the Rabi frequency tt gr has been rescaled by 
1 / VN to correct for the y/N enhancement due to the Ry- 
dberg blockade effect. From the figure it is clear that the 
Gaussian 7r-pulscs fail for Lo > 0.7 /j,m (corresponding to 
more than three atoms) where fl gr gets comparable to 
the interaction strength. On the contrary the optimized 
pulses yield infidelity e < 10 -2 up to Lq ps lA/im, thus 
allowing to almost double the system size. So by opti- 
mization we can outperform the blockade radius associ- 
ated with the guess pulse. This effect exceeds what one 
can expect by just stretching tt gr over the whole opera- 
tion time and thus lowering its bandwidth. The presented 
pulses have only real Rabi frequencies. Allowing also for 
complex values (via phase modulation) roughly results 
in a 15% relative improvement of the state preparation 
error. 

As stated in the introduction, our final goal is to pre- 
pare an ensemble of atoms with (uniformly distributed) 
random positions in 3D space moving randomly with 
thermally distributed velocities in the W-state. We thus 
now consider a 3D cloud of N frozen atoms at random 
positions. Differently from the previous case, the inter- 
actions between atoms are random due to the distance- 
dependent interaction terms and this means that we have 
to produce very robust pulses "on average" at the cost of 
some fidelity loss for each given sample of atoms. To this 
aim, we extract the relevant features of the optimized 
pulses e.g. Fig. [3] (middle) and optimize only few param- 
eters like the height and the width of the two pulses. The 
resulting optimal pulses are reported in the right panel of 
Fig. [3] improvement here is mainly due to the reduction 
of the bandwidth. The \g) <-> \r) transition is performed 
via a long fiat pulse followed by a fast kick in the \r) <H> |e) 
transition. Finally, we compare the performance of the 
two schemes: the infidelity e obtained by W-state prcpa- 
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FIG. 4. 3D cloud of 10 randomly distributed atoms: infidelity 
e as a function of the cloud's diameter Lo plotted for Gaus- 
sian 7r-pulses (circles) and an optimized sequence (squares). 
The inset shows that N has almost no influence on e. The 
curves in the inset correspond to Lo — 0.83. . . 1.46 /im and 
the optimized pulse. 



ration with the optimized pulse and with the guess Gaus- 
sian 7r-pulses averaged over different system samples (i.e. 
different instances of the random positions). Fig.|4]shows 
the infidelity e for clouds of ten atoms as a function of 
the maximum cloud extension Lq while in the inset of 
Fig. H e is shown as a function of the number of atoms 
N = 7, 8, ... 11 in the sample. In both cases the errorbars 
correspond to the statistical noise from 8 realizations of 
random atoms positions within a cloud of fixed maximum 
diameter Lq. The optimized pulses clearly result in im- 
proved infidelities and most importantly with almost no 
dependency on the number of atoms. This result, as we 
shall show below, is corroborated also by our analytical 
estimate of accuracy of the optimal state preparation. 
An additional complexity aspect for the experimental re- 
alization of the proposed protocol arises from the fact 
that the number N of atoms in the ensemble depends 
on the density and the excitation volume and can only 
be determined with Poisson precision (±y/N) changing 
a 7r- pulse roughly as cos(^Vl ± N) pa 1 — This 
very weak dependence on the atom number N of the fi- 
nal fidelity paves the way to a successful experimental 
realization of the W-state at room temperature. 

We can provide an analytical estimate of the state 
preparation error under the hypothesis that: a) the major 
deviation of the prepared state from the desired W-state 
is given by the population in the two-excitations sector, 
b) each Rydberg atom is approximately performing an 
independent Rabi oscillation between \g) and \r) and c) 
the influence on an atom by another nearby excited one 
is to detune the single atom dynamics. The detuning will 
just be the van-der-Waals interaction between two Ryd- 
berg excited atoms V = Ce(L~ e , with dij the distance 
between atom i and atom j and thus the resulting Rabi 
oscillation for this single atom is then only going up to 
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FIG. 5. 3D cloud of N = 10 randomly distributed atoms: 
Errors in the state preparation. Pi calculated by MPS simu- 
lation is compared to the theory curve from the P(r) model. 



having a fraction of 
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in the excited 



W gr +V 2 ~ 1 l+c 2 dYj 

(Rydberg) state (c = fl gr /Ce). This double excitation 
part is missing in the population Pj representing the first 
excitation sector. Summing this up over all neighbors 
gives 
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In the thermodynamical limit an analogous continuous 
expression can be obtained by considering a homogeneous 
spherical cloud of radius R, introducing P(r) as the single 
atom excitation. This yields 



P(r) 



1 

NV 



p 2 sin OdpdOi 



1 + c 2 (r 2 + p 2 — 2rp sin cos < 



(5) 



with r the radial position (distance from the center of 
mass of the cloud) of atom i in the cloud and integrating 
over the sphere's volume. Fig. [3] shows the comparison 
between Eq. (O and the Pi obtained from a numerical 
simulation of the dynamics of ten atoms. The good corre- 
spondence of this high density limit theory and numerical 
results of the few body simulation supports the previous 
finding of the very weak dependence of the results on the 
number of atoms TV, that is, considering only ten atoms 
already is enough to get the major features of the state 
preparation. 

If we now include in our theoretical description the fact 
that the atoms are moving with a thermal distribution, 
we shall consider the lasers' Doppler shifts depending on 
the velocity of the atoms. Furthermore we shall analyze 
the effects of a time-dependent interaction V(t) and the 
fact that the atoms might move out of the laser beam 
during state preparation, as at room temperature an av- 
erage atom moves about 0.5 /im in To time. For this time 
scale we can achieve good blockade within a sphere of ra- 
dius 0.5 — 0.55 ^m (corresponding to Lq = 1 — 1.1/im). As 
the pulse is robust with respect to random positions we 
expect that the thermal motion will not affect drastically 
the results. 

Finally, the Doppler effect will cause a widening of the 
excitation populations with a Lorentzian shape Pf oc 
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k /ft ) depending on the velocity ivii of atom 
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FIG. 6. Directionality p as a function of temperature ■& for 
parallel lasers (red) and anti parallel lasers (black). Squares 
are values from state preparation with optimized pulses, cir- 
cles from Gaussian 7r-pulses. The inset shows the emission 
cone on resonance for $ = 220° C and anti parallel lasers. 



(4) shift enters as an additional term (k% ± fe^y \r)(r\ into 



in the direction of the incoming lasers. The Doppler 



Hioc,i (we neglect the Doppler shift of the third laser since 
f2 re 3> k re Vi). The wave numbers ki of lasers 1 (|<?) <-> \i) 
transition) and 2 -s-> |r) transition) are summed up or 
subtracted depending on whether the lasers are parallel 
or anti parallel. We consider a sphere of radius 0.53 /im 
and an atomic density and velocity distribution corre- 
sponding to temperatures of 200 — 260 °C, that is approx- 
imately N = 200 — 1200 atoms and a velocity distribu- 
tion characterized by a Gaussian width of 213 — 226 m/s. 
We calculate the approximate W-state \W) = J2i a i\ e i) 
(where the cv^ contain the Doppler widening depending on 
the directions of the lasers) obtained by our state prepa- 
ration pulses by simulating the evolution by fl gr for N 
atoms assuming perfect blockade, thus keeping track only 
of the ground state and the N singly excited states and 
constraining the system evolution to the ground state and 
the first excitation sector. To model the decay of \ W) we 
follow the full exponential kernel description [li|ll|3l| 
and in addition we include also the motion of the parti- 
cles. 

Fig. [S] shows the resulting directionality p (the prob- 
ability of photon emission in the forward direction) as 
a function of the temperature, each value averaged over 
240 random realizations of a uniform position distribu- 
tion and Gaussian velocity distribution after a decay time 
of 100 ns. The inset shows the emission cone on reso- 
nance for i9 = 220 °C. The forward cone is defined to 
cover about 3% of the solid angle corresponding to a 
maximum deviation from the forward direction of 0.3 rad 
(dashed lines) . Values for p are black for parallel and red 
for anti parallel lasers on the \g) <-> \i) and \i) <-> \r) tran- 
sitions. Clearly, the highest directionality is obtained for 
anti parallel lasers. Squares correspond to state prepara- 
tion with the optimized pulse. Circles mean state prepa- 
ration with Gaussian 7r-pulses - in this case we have to 
decrease the cloud radius and thus deal only with about 
half as many atoms participating in the W-state leading 



5 



to a reduction of the directionality. 

At higher temperatures samples are denser and thus 
more atoms can be excited, resulting in higher direction- 
ality. However, the collision rate is enhanced thus co- 
herence times are decreased. Therefore the temperature 
has to be chosen properly. The best choice in the exper- 
imental settings is assumed to be in the range of 220° C, 
corresponding to an overall 70 % increase in directionality 
compared to the guess pulse. 

Conclusions - We have analyzed the possibility of 
constructing a single photon source using a thermal vapor 
of Rubidium 87. The results, however, are also qualita- 
tively valid for other atoms like Cesium. By CRAB as- 
sisted pulse shaping using a Matrix Product State code 
we showed how to increase the number of atoms in a 
W-state by a factor of two compared to the initial guess 
which corresponds roughly to an increase of 70 % in di- 
rectionality of the emitted photon. Choosing the exciting 
lasers as anti parallel gives an additional crucial increase. 
We plan to experimentally test these predictions also to 
quantify the effect of the collisions and laser focus not 
included in the present description and to see how they 
compare to experiments with ultracold atoms 
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